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ERRATA 
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1)  page  15,  line  9  from  bottom,  read  "dt"  for  "dy". 

2)  page  18,  line  4  from  bottom,  read  "(23a)"  for  "(23aa)". 

;5)  page  J>k ,   line  8  from  bottom,  read  "any  plane  x  =  x  with 
|x^-  x|  =  o(l/e)  at  a  time  t  =  t" , 

k)     page  35,  line  J,   read  "lim  ev  <  0". 

5)  page  36,  line  6  from  bottom,  "see  Lemma  Al,  in  Appendix  3". 

6)  page  38,  line  7,  read  "positive  value  of  u"  for  "positive 
value  u" . 

7)  page  39'   In  the  second  figure  the  axes  are  x  and  u. 

8)  page  45.   In  the  figure,  the  curve  cuts  axis  orthogonally 
at  X  =  X,  . 

9)  page  55,  line  5  from  bottom,  read  "r^  =  [(u-l)^+  v^L  ". 


NYO-2885 
PHYSICS  AND 
MATHEMATICS 

ABSTRACT 

In  this  report  we  present  a  method  for  studying  the  internal 
sti*ucture  of  a  magneto -hydrodynamic  shock  in  which,  in  contrast 
to  earlier  theories,  there  are  no  collisions  among  the  plasma 
particles.   The  theory  yields  us  two  distinct  situations.   In  the 
first  case,  there  is  no  shock.  The  ions  and  electrons  pass 
through  a  region  of  changing  magnetic  and  electric  field  and 
eventually  they  return  to  the  state  from  which  they  started.   In 
the  second  case,  there  is  a  definite  change  of  state.  Within  the 
approximation  we  obtain  a  qualitative  pictxAre  of  the  shock.  The 
magnetic  field  rises  and  then  oscillates  with  a  characteristic 
wave  length  about  a  mean  value  greater  than  its  initial  value. 
The  distribution  function  of  the  ions  is  greatly  distorted.  The 
approximation  is  based  on  an  asymptotic  development  in  the  ion  to 
electron  mass  ratio.   It  is  valid  over  distances  that  are  comparable 
to  or  even  large  compared  to  the  wave  length  of  the  oscillation  but 
small  compared  to  the  ion  Larmor  radius.  Thus  it  is  quite  possible 
that  an  improvement  in  the  range  of  validity  of  the  approximation 
would  show  that  the  oscillation  is  damped  out  (e.g.,  on  a  length 
scale  comparable  to  the  ion  Larmor  radius)  and  that  a  classical 
shock  exists.  However,  the  change  of  state  shows  that  an  irrevers- 
ible process  has  been  demonstrated  to  be  possible  without  collisions, 
and  this  already  marks  the  difference  between  existence  of  shocks 
and  nonexistence. 

It  should  be  noted  that  the  fields  and  flows  are  self-consistent 
within  the  approximation. 
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NYO-2885 
MAGNETO  HYDRODYNAMIC  SHOCK  STRUCTURE  WITHOUT  COLLISIONS 

1.   General  Introduction. 

To  determine  the  Internal  structure  of  a  magneto-hydrodynamlc 
shock  one  wants  to  establish  the  existence  of  a  steady  state 
solution  of  the  one -dimensional  equations  of  motion  for  two 
fluids  and  Maxwell's  equations,  with  the  property  that  It  leads 
from  a  given  constant  state  to  a  constant  state  at  higher  pres- 
sure.  If  the  second  state  Is  the  same  as  the  first  we  call  the 
solution  a  pulse.   This  problem  has  been  under  study  for  some 
time.   Early  theories,  see  [11],  [12],  [14],  [17],  all  Involved 
In  some  sense,  collisions,  either  through  viscosity,  thermal  con- 
ductivity or  electric  conductivity.   In  those  cases  where  a 
transition  stage  leads  from  one  constant  state  to  another  the 
length  of  this  state  is  rather  long  —  of  the  order  of  a  mean  free 
path.  This  distance  may  actually,  in  certain  cases,  compare  with 
the  dimensions  of  the  apparatus.   For  a  useful  shock  theory  one 
needs  a  formulation  that  does  not  Involve  collisions.  This  is  a 
radical  departure  from  the  usual  procedure  and  it  was  first  proposed 
by  H.  Grad,  see  [7]  and  [9],  that  such  a  theory  was  possible.   Also 
such  a  shock  theory  can.  In  contrast  to  a  collision  theory,  lead  to 
a  difference  in  the  temperature  of  ions  and  electrons. 

We  begin  by  assuming  that  there  are  no  collisions  and  do  not 
treat  the  gas  as  a  continuous  fluid.   We  seek  a  flow  and  field  that 
are  self-consistent;  this  means  that  the  influence  of  ions  and 
electrons  on  electric  and  magnetic  fields  and  vice  versa  is  taken 
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into  account  .   Properly  speaking,  we  seek  a  time  independent 
one -dimensional  solution  of  two  Boltzmann  equations  (one  for  ions 
and  one  for  electrons)  and  of  the  Maxwell  equations  for  the 
electromagnetic  field. 

If  x'  is  the  space  variable,  the  conditions  at  x'  =  -»  are 
that  the  fields  are  given  constants  and  the  velocity  distributions 
of  the  ions  and  electrons  are  also  given.   We  consider  the  follow- 


— » 

E 


X  ' 


direction  of  mean 
flow 


Ing  situation:  the  electric 
field  is  in  the  y' -direction  and 
the  magnetic  field  is  in  the  z'- 
direction.  The  distribution  func- 


tions f  for  the  ions  and  f  for 


the  electrons  are  prescribed  to 
be  compatible  with  a  constant 
state.   As  we  shall  show,  we  can  work  with  distribution  functions 
at  -00,  f±   which  depend  only  on  u',  v',  the  velocity  components  in 
the  X  and  y  directions.  We  take  the  distribution  function  for  the 
ions  at  -»  as  a  function  of  (U-u' )  +  v'  where  U  is  the  mean  velo- 
city at  -00.   Since  the  flow  is  assumed  to  be  steady,  U  is  also  the 
shock  speed. 

In  order  to  simplify  the  problem  we  take  f  =  0  for 
(U-u')^+  v'^  >  R'^.   We  refer  to  R'  as  the  cut-off  radius.   Depend- 
ing on  how  we  choose  the  cut-off  radius  we  can  observe  different 
effects. 


A  totally  different  approach  to  this  problem  involving  "turbulence" 
in  the  magnetic  field  is  given  in  [6].   A  fluid  approach  involving 
strong  shocks  is  given  in  [4]. 
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The  fields  at  -oo  are  restricted  to  be  those  for  which  the 
Alfven-Mach  nimber  1  <  U/a  <  v  <  2  where  a  is  Alfven  speed  and 
V  is  some  number  depending  on  the  distribution  at  -co.   Its  value 
must  be  determined  in  each  case  by  numerical  computations. 

The  basis  for  approximation  is  that  e  =  m  /m  ,  the  ratio  of 
electron  mass  to  ion  mass,  is  a  small  number.  We  expand  in 
powers  of  e  in  all  equations  keeping  an  appropriate  scale  length 
fixed.  This  length  is  \/ 1  I  where  I  is  the  distance  a  particle 
moves  forward  in  a  complete  change  of  velocity  phase,  that  is,  a 
gyro-period,  (see  Appendix  l). 

One  of  the  keys  to  finding  a  solution  by  an  expansion  proce- 
dure is  to  take  E„, ,  the  charge  separation  field,  of  the  form 
E(x',€)/e  where  E  is  finite  as  e  — p>  0.   In  other  words,  there  is 
a  large  charge  separation  field  due  to  the  disparity  in  the  masses 
of  the  ions  and  electrons.   Another  vital  point  is  to  note  that 
for  €  ?^  0,  no  matter  how  small,  all  ions  travel  from  -»  to  -h» 
provided  H  ,  >  0.   It  will  be  shown  that  as  c  — >0  this  requires 
that  special  care  must  be  used  in  determining  the  asjmiptotic 
behavior  of  the  ions,  since  there  are  no  permanently  trapped  ions 
and  the  behavior  to  lowest  order  in  e  seems  at  first  glance  to 
imply  that  there  are.   In  computing,  we  have  to  make  the  assump- 
tion that  there  is  a  point  x'  =  x^  such  that  for  x'  >  x^  the 
charge  separation  field  dominates  and  for  x'  <  xj[.  the  magnetic 
field  dominates. 

The  results  which  are  proved  in  the  body  of  this  report,  for 
lowest  order  in  e,  are  as  follows:  There  is  a  number  R '  <  U  such 
that 
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I.   If  R'  <  RS  all  ions  go  from  -»  to  +<»  without  t\irning 
around.  The  magnetic  field  has  a  pulse  but  all  quantities 
return  at  +»  to  the  state  they  were  in  at  -»,  see  Figure  1. 
This  solution  confirms  to  lowest  order  in  e  certain  special 
results,  the  well-known  limiting  solution  for  zero  pressxire 
(R'  =  0),  see  Section  4  and  also  [1],  and  Mittleman's  unpublished 
results.  There  is  no  shock  and  the  flow  is  completely  reversible. 
Simple  two-fluid  theories  (both  adiabatic)  without  standard  dis- 
sipative  mechanisms  lead  to  similar  flo\i?s  (not  shocks),  see  [1], 
[7]>    [8],  [12].   Refined  two-fluid  theories  have  so  far  also  not 
led  to  shocks.   According  to  Liist '  s  approximation,  a  two-fluid 
theory  with  non-symmetric  pressures  does  not  yield  a  shock  solu- 
tion, see  also  [15]. 


^x' 


Hz' 


X' 


Figure  1 


II. 


R  '  <  R'  <  U.   Some  ions  bounce  a  few  times.   In  doing 


so  they  increase  the  density  in  front  of  the  magnetic  pulse.  As 
a  result  the  symmetry  about,  say  x'  =0,  disappears  and  the  solu- 
tion behaves  for  x'  <  0  like  thg  first  part  of  the  single  pulse 
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but  for  x'  >  0  it  oscillates  without  damping.   See  Figure  2, 

H 


Figure  2. 


2.   Equations  and  Bovdary  Conditions. 

The  equations  of.  our  system  are  Maxwell's  equations  for  a 
steady  state  electromagnetic  field  and  two  Boltzmann  equations 
without  collision  terms.  The  former  are 


(1) 


/ 


< 


curl  H  =  J 
curl  E  =  0 
dlv  ixH  =  0 
div  K'E  =  q 


where  H  is  magnetic  field,  E  is  electric  field,  J  is  current, 

it 
q  is  net  charge  per  unit  volume,  K'  and  ij,  are  vacuum  constants  . 

The  two  Boltzmann  equations  reduce  to 


(2) 


""i  ^  •  ^x'^+  t  e(E  +  u  X  M.H)  •  7^,  fj.  = 


since  there  are  no  collisions.  Here  m^  is  the  mass  (of  the  ion 


The  units  in  the  above  equations  are  m.k.s, 


9.  - 


or  electron),  u  =  (u',v',w')  is  the  velocity  vector.  The  space 
variables  are  x',y',z';  e  is  the  electron  charge  and  f_|^(u' ,v' ,w' ,x' ) 
is  the  distribution  function  (for  ions  or  electrons),  i.e.,  the 
niJinber  of  ions  in  an  element  of  phase  space  per  unit  volume. 

The  current  and  charge  in  Maxwell's  equations  (l)  are  related 
to  the  distribution  functions  f '  by 
r  A 


(3)    < 


J  =  e  ///(f|  -  fMu  du'dv'dw' 


q  =  e  ///(f|  -  f^du'dv'dw' 


where  the  domain  of  integration  is  the  whole  space. 

We  assume  that  all  quantities  depend  on  one  space  variable 

only,  say  x".   As  a  result  we  see  from  (1)  that  H  , ,E  , ,E  ,  are 

X   y   z 

all  constant.  We  take  H  ,  =  E  ,  =  0.   We  prescribe  H  .  =0  at  -co 

X     z  y 

and  we  look  for  a  solution  with  H   =0.  The  two  varying  fields 

y' 

are  E  ,  and  H  ,  and  we  prescribe 


(^) 


L 


^x'  =  0  >   Hz'  =  K^ 

f»  =  F'(u',v',w')    for  X'  = 

X       X 


The  prescribed  distribution  functions  must  be  consistent  with  a 
constant  state. 

The  problem  of  finding  a  strict  shock  solution  is  to  find  a 
solution  of  (l)-(4)  that  tends  to  a  constant  state  as  x'  — > +« 
and  which  differs  from  the  state  at  -co.   No  such  solutions  have 
been  found.   In  fact  no  exact  solutions  to  (l)-(4)  except  the 
solutions  with  zero  pressure  have  been  determined.   Attempts 
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to  find  solutions  nimerlcally  by  forward  differences  in  x'  have 
been  unsuccessful  because  the  characteristics  of  (2)  in  the 
u' ,v' ,w' ,x' -space  may  intersect  the  planes  x'  =  constant  more 
than  once.  The  present  method  for  studying  solutions  is  to 
perform  asymptotic  expansions  in  a  suitable  parameter  with  appro- 
priately normalized  variables. 

But  first  we  introduce  some  notation  and  derive  certain 
conditions  for  x'  — > -»  which  are  useful  in  normalizing  the 
variables.  We  introduce  the  nianber  densities 

(4a)      n.  =  //  f|  du'dv'dw' 

..,1':  , 

and  the  mean  velocities 

u!  =  //  u'f'  du'dv'dw'///  fl  du'dv'dw' 

?•  =  //  v'f'  du'dv'dw'///  f'  du'dv'dw'  . 

Integrating  the  Boltzmann  equations  (2)  we  find 

(4c)      '^+^-1  ^^^  ^J^l     ^^®  constant. 

Since  the  fields  tend  to  constants  as  x'  — > -oo  we  see  from 
(1)  that  q  — »  0  and  J  — $>0.   From  (3) 

J  =  e(n  u'  -  n  u' )  and  q  =  e(n,-  n  ), 

and  hence  n  -  n_  — $>  0  as  x  — >«,  and  u'  -  ul   — >  0.  Thus  setting 
n  =  n*  and  u'  =  U  at  -«,  we  have 
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(5)  n  =  n_        for  x  =  -co 

(6)  u '  =  u^  =  U       for  X  =  -00 

and 

(7)  n_^u|  =  n_u^  =  n"u     for  x  =  - 


Note  that  U,  the  mean  value  of  the  speed  of  the  Ions  and  the 
electrons  as  x'  — >  -oo,  is  the  speed  of  the  shock  viewed  from 
a  frame  of  reference  in  which  the  mean  speed  of  the  gas  is  zero 
and  the  shock  is  moving. 


Prom  (2)  we  have  as  x'  — >  -co,  since  f'  become  independent 


of  X', 


(E  +  u  x  H)  .  Vy'^i  =  ° 


and  hence  using  (4)  and  (6)  we  find  that 
(7a)      Ey,  =  nUH^,. 

We  now  introduce  an  intetrated  distribution  function 

(7b)      f^.  =  I  f;  clw'. 


-00 


Then  equation  (2)  holds  for  f  instead  of  f !.  The  next  step  is  to 
normalize  the  variables.  . — :;;;;-. 

We  want  to  expand  with  respect  to  the  square  root  of  the  mass 

ratio,  and  we  find  that  the  appropriate  length  to  keep  fixed  is 
given  by 


(8)       X  =  e 


CO 

n  n 

-i-  X' 
m 
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The  other  suitable  dlmensionless  variables  are 


(9) 


with 


u 

= 

u'/U 

V 

= 

v'/U 

H 

= 

«./«:• 

v 

E 

e 

CE^,/Ey, 

ym_/in_^ 

Note  that  with  this  choice  of  variables,  E  ,  — >  ~,  as  e  — >  0. 
Using  (7)  and  (3),  Maxwell's  equations  (1)  reduce  to 


(10) 


K 


dE 
dx 


1_ 
u, 


1_ 
u 


2  / ^ 

where     K  =  ^—ik  ^     and  A  is  the  Alfven-Mach  number,   A  =\/ra  n     U  /p-lCi 
Ag2  ^^  y    +  +00  z 

c  =  l/  VK'ii  is  the  speed  of  light,   v     is  the  mean  velocity  in  the 

y-direction  after  normalization,   that  is; 


(11) 


-1 


V.  =  //  f ,v  du  dv/J/  f,du  dv  =  n"  //  f  v  du  dv, 

+    ■'■'+  ±  +i- 


and  u  is  the  normalized  mean  velocity  in  the  x-direction,  that  isj 


We  note  that  K'  is  about  10"^.  Hence  if  u/c  is  sufficiently  small 
compared  to  €^  so  that  K  is  small  compared  to  1,  we  may  take  u+=  u_. 
In  this  case  the  total  charge  is  zero.  This  approximation  which  is 
frequently  made  will  be  valid  unless  dE/dx  turns  out  to  be  very 
large.  However,  it  is  not  necessary  to  make  this  assumption,  and  we 
do  not. 
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(12) 


-1 


u  =  //  f_^u  du  dv  /  //  f  du  dv  =  n^"^  //  f  u  du  dv 


For  a  suitable  choice  of  the  distribution  functions  at  -«  we  have 
V  =  0  at  -00.   On  the  other  hand  by  (4a),  (4b)  and  (4c),  (6)  and 
(7b), 

(12a)     //  f_^u  du  dv  =  1. 

With  the  change  of  variables  (7b),  (8),  (9),  the  two 

Boltzmann  equations  (2)  reduce  to 

r 


(13) 


f^f  f^f  r*f 

Au  ^  +  (E  +  evH)  ^  +  €(1-UH)  ^  =  0 


;^f 


Sf 


sf 


Au 


^  -  ?  (|  ^  -«)  ^  -  ?  (1-^)  ^  =  °- 


Vx 


From  the  boundary  conditions  (4)  at  -co,  we  have  with  (9) 


(14) 


H-1  =  E  =  0,   atx  =  -ooj 


(15^)     f^  =  C  , 
(15_)     f_  =  f!  , 


at  X  =  -oo; 


at  X  =  -00. 


We  restrict  f  to  satisfy 

f  f *  =  P((l-u)^+  v^)  ^   0  for  (l-u)^+  v^  <  R^ 


(16) 


=  0 


for  (l-u)^+  v^  >  R^, 


The  cut-off  parameter  R  =  R'/U  will  be  specified  later.  We  Impose 
here  only  the  condition 


(I6a) 


R  <  1, 


14  - 


This  means  that  for  x  — >  -«  the  ions  all  have  positive  velocity 
in  the  x-direction  and  a  speed  less  than  2.  The  function  f** 
may  be  any  monotone  function  of  speed  at  -e». 

The  problem  is  to  solve  (10),  (11),  (12),  (15)  under  the 
conditions  (l4),  (15.),  (l6)  and  (l6a).   We  require  that  E,H 
are  bounded  and  have  bounded  derivatives  for  all  e. 

We  have  in  addition  from  (11),  (12),  (6),  (l6) 


(l6b)     u  -  1  =  u_-  1  =  V  =  v_  =  0   at  X  =  -00 


From  (13)  we  see  that  the  function  f  is  constant  along  the 
electron  particle  paths: 

r 


(17)    \ 


du  _   du  _    1  /E    „x 


|v_dv_^(,.^, 


where  dy  =  dx/u  is  the  differential  of  the  time  parameter  on  a 
particle  path. 

On  the  other  hand,  the  function  f  is  constant  along  the  ion 
particle  paths: 


(18)    { 


dU      ,,   dU      1   /ir,   .      tt\ 


V, 


dv 
dT 


dv 


=  uS=|(l-.H). 


The  method  for  finding  an  approximate  solution  is  the  follow- 
ing. First  we  find  asymptotic  expansions  in  e  for  (ll)_  ajJid  (12)_ 
for  the  electrons.   It  is  then  possible  to  find  integrals  of 
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equations  (10)  and  (15)  and  we  are  left  with  a  Boltzmann  equation 
for  the  ions  with  coefficients  which  depend  on  moments  of  the  ion 
distribution  function,  (Section  5).  We  then  look  separately  at 
the  various  characteristics  (l8)  for  the  ion  Boltzmann  equation 
and  their  dependence  on  x  for  € — >  0,  (Sections  5  and  7)  and  thereby 
draw  conclusions  about  the  behavior  of  E  and  H  (Sections  6  and  8). 


3.   Asymptotic  Behavior  of  the  Mean  Velocities  of  the  Electrons 
and  Resulting  Momentum  Integrals. 

The  second  equation  of  (13)  for  the  distribution  function 
for  the  electrons  will  now  be  solved  for  E  and  H  by  an  asymptotic 
series  in  e,  and  from  this  expansion,  the  average  values  of  u_ , v_ 
which  occur  in  (10)  may  be  determined.   We  have  then  a  simplified 
problem  to  solve. 

However,  it  is  not  possible  to  find  directly  from  (13)  an 
expansion  of  f_  in  a  power  series  in  e .   A  preliminary  transforma- 
tion is  necessary  and  this  transformation  is  suggested  by  an 
approximation,  very  similar  to  the  "guiding  center"  approximation 
but  differing  from  it  because  of  the  large  charge  separation  field. 

We  note  that  if  we  set  v  =  -  E/eH,  u  =  1/h  +  ^  (^)  in  (1?) 
the  equations  are  satisfied  to  the  highest  order  in  €.   This 
suggests  Introducing  as  new  independent  variables 


< 


^'  =  -  H  ^  ^  (M)  ^  " 
"  ^  d3F  ^H  ^ 


(19) 

■n'  =  +  eAh  +  V. 
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With  these  nei>r  variables  it  is  possible  after  some  calcula- 
tions, which  are  given  in  Appendix  2,  to  find  an  asymptotic 
expansion* for  f  .   Substituting  this  expansion  in  (11)  and  (12), 


one  finds  to  order  e 
u  =  3 


(20) 


ev 


"  +  d5^  ^H  ^ 


=  -  e/h. 


L 


We  might  interpret  this  by  saying  that  in  the  mean  the 
electrons  behave  like  a  guiding  center  fluid.  There  exists  an 
adiabatic  invariant  for  this  motion  and  it  is  therefore  possible 
to  show  that  the  electrons  are  heated  adiabatically.   (This 
adiabatic  invariant  is  not,  however,  the  magnetic  moment.) 

Substituting  (20)  in  (10)  we  find,  along  with  (13),  that  the 
problem  is  reduced  to  solving 

r 


^    (H  +  I  n^)  =  TlH  -  ^r^ 


asE" 


Aev, 
u, 


(21'     l^(i^¥'  =  -^  + 


H 


\ 


^^^+  (-  J  +  cv)H^+e(l-uH)^=0 


where  u     and  v     are  given  by  (11 )  and   (12)  and  tj  =  -AE/h.     The 
boundary  conditions  reduce  to 


H  -  1  =  T)  =  0       at  X  =  -co 


(22) 


f.   =   C 


For  a  general  theory  of  such  asymptotic  expansion  see  a  forth- 
coming report  by  C.  S.  Gardner.         , 
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Fields  E  =  -  tjH/a  and  H  which  satisfy  (21)  and  (22)  will  satisfy 
(10-15  )  to  terms  of  order  e  . 


Energy  and  momentimi  integrals  of  the  equations. 

Two  Integrals  of  equations  (21)  can  be  found  corresponding 
to  the  conservation  of  momentvun  flux  In  the  x  and  y  directions. 
Multiplying  the  last  equation  of  (21)  by  u  and  Integrating  by 
parts  over  the  (u,v) -space  we  find 

A  ^  //  f^u^du  dv  +  fj   Q  -   €v)Hf_^du  dv  =  0 

or  using  the  first  two  equations  of  (21),  with  (11),  (12)  and 
(12a): 

2 

(23a)     I(x)  +  5_p-  (1  -  Kti^)  =  constant 

2A 
where 

(23aa)    I(x)  =  //  f_^u^du  dv. 

The  constant  Is  foimd  using  (22),  to  be  I(-oo)  +  l/2A  . 

Similarly,  one  finds  by  multiplying  the  Boltzmann  equation 
by  V, 

(23b)     AJ(x)  +  €Ti(l  +  ^)  =  0 

where  J(x)  =  //f  uv  du  dv. 

It  can  be  shown  that  (23aa)  and  (23b)  may  be  solved  for  t^  and 
H  in  terms  of  I  and  J.  The  functions  I  and  J  are  simply  related 
to  the  stress  tensor. 

An  additional  Integral  relation  involving  the  electric 
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potential  6  =  -!   Edx  proves  useful  In  the  lowest  order  case.  From 
the  first  equation  of  (21)  we  note,  "by  integrating  from  -»  to  x, 
that 


(23c) 


AJ?5  =  -A  I  E  dx  =  I  TiH  dx  =  -1  +  H  +  -%-  +  €A 


2  .-7, 


2 


U. 


.2       r^v 


=   constant  +H  +  ^+€A    -idx  for  any  x^,.. 


4.   The  Simplified  Problem 

The  original  boundary  value  problem  has  been  reduced  to 
solving  the  Boltzmann  equation  for  the  ions,  (the  last  equation 
of  (21)): 

(24)      Au  ^  +  (-  J  +  ev)H  ^  +  e(l-uH)  ^  =  0 

where  t)  and  H  are  given  in  terros  of  the  moments  l(x)  =  //  f  u  du  dv 
and  J(x)  =  //  f  uv  du  dv  through  (25a)  and  (25b)  with  one  condition 
prescribed  at  -»,  namely,  f  =  f°°  where  f^  is  given  by  (l6).  We 
also  have  (25c)  which  turns  out  to  be  easier  to  use  than  (25b). 
We  would  like  to  know  when  there  exists  a  solution  to  this 
problem  aside  from  the  trivial  solution  f_^(u,v,x)  =  f°^(u,v).  First 
of  all  one  can  show,  (see  Calkin  [5]  for  the  general  proof  )  that 
there  exists  a  solution  near  x  =  -w  where  f  increases  exponentially 
in  X.   Then  by  (25a)  and  (25b),  tj  and  H  change  exponentially  also. 
We  want  to  know  if  this  solution  can  be  continued  in  x. 
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In  accordance  with  the  general  procedure  so  far  we  would  like 
next  to  find  an  asymptotic  expansion  of  (24),  (23a),  {^yo) ,   with 
the  "boundary  condition  (l6),  to  lowest  order  in  e.  However,  we  run 
into  two  difficulties  if  we  look  only  at  the  lowest  order  tex*ms. 
First  of  all,  it  turns  out  that  we  cannot  solve  to  lowest  order 
and  satisfy  the  boundary  condition  at  -oo,  (l6),  in  a  straight- 
forward fashion.   Secondly,  as  we  shall  see,  to  lowest  order  in  €, 
there  are  closed  characteristics  and     " 
since  they  plainly  never  reach  -co  it 


seems  that  we  could  prescribe  f       ^ 
arbitrarily  on  them.   In  this  case 


^ 


the  solution  does  not  appear  to  be 
* 


Lowest  order  characteristics 
in  x,u  plane. 


unique 

However,  these  difficulties  can  be  overcome  by  using  the 

second  order  terms  in  the  expansion  and  some  general  arg\jments 

about  the  behavior  of  the  characteristics.  The  correct  procedure 

is  to  use  the  lowest  orde v*  term  in  the  expansions  in  e  to  decide 

on  what  curves  f  ir^  cone  pant  and  to  use  more  refined  arguments 

to  decide  what  value  of  f  to  associate  with  these  curves.  The 

major  part  of  the  work  is  concerned  with  the  latter  problem.   Once 

we  know  to  lowest  order  in  €  the  curves  on  which  f ,  is  constant 

+ 

and  qualitatively  what  the  constant  is,  we  can  draw  conclusions 
from  (23a)  and  (23c)  about  the  behavior  of  H  and  r\. 

In  the  next  section  (5)  we  shall  discuss  the  very  simple  case 
of  zero  pressure  which  gives  us  some  insight  into  some  general 


Sec  [2]. 


-  SO  - 


properties  of  the  solution  of  (24).  Then  In  Section  6,  we  retiim 
to  the  problem  of  finding  the  lowest  order  characteristics  and 
attaching  values  of  f  to  them. 

In  the  rest  of  this  section,  we  shall  draw  some  conclusions 
about  characteristics  in  general.  Now  we  consider  general  pro- 
perties and  characteristics.  We  use  Interchangeably  the  words 
"particle  path"  and  "characteristic". 

We  recollect  that  f  is  constant  on  the  characteristics  which 

are  given  by  (l8).  Now,  if  for  all 

X  >  X  we  have  u  >  0  on  the  particle 

paths  issuing  from  x  =  x  ,  f  is 

determined  for  x  >  x  by  its  values 

o 

on  X  =  X  .   Now  suppose  that  some 

particle  paths  loop,  i.e.,  u  cha:;\ges 

sign  at  X  =  x, ,  say,  and  a  particle  path  comes  back  to  x  =  x  . 

Then,  on  this  particle  path 
a  certain  value  of  f  is 
carried  along,  see  figure. 


and  plainly  f  cannot  be 
arbitrarily  prescribed  on 


x=x 


X=Xn 


X  =  x  . 
o 


c  1 

However,  in  our  problem  under  certain  simple  and  reasonable 

assumptions  about  the  field  it  will  be  shown  that  no  particle 
path  will  return  to  x  =  x  if  x  is  sufficiently  large  and  nega- 
tive. Therefore,  we  conclude  that  f"*"  is  determined  everywhere 
by  its  values  for  large  negative  x. 
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The  necessary  assimiptlons  are  that  \l6\   =   |-/  E  dx|  is 
bounded  and  that  H  is  bounded  and  bounded  away  from  zero.  These 
conditions  hold  for  all  approximate  solutions  that  are  foiond 
here.  The  above  statement  about  the  behavior  of  the  ions  is  a 
special  case  of  a  more  general  lemma. 
Lemma  1%     Every  characteristic  of  (24)(given  by  (l8))  extends  from 

X  =  -00  to  X  =  +00. 

This  lemma  is  equivalent  to  the  statement  that  there  are  no 
trapped  ions  and  will  be  proved  later. 

A  second  lemma  will  be  needed  which  tells  us  how  v  increases 
with  time,   A  characteristic  which  loops  and  crosses  some  plane 
X  =  constant  several  times  will  have  different  values  of  u  and  v 
associated  with  each  time  it  crosses  such  a  plane. 
Lemma  2;  On  a  characteristic,  the  values  of  v  corresponding  to 
the  same  value  of  x  increase  with  time. 

Proof;  Along  any  characteristic  we  may  introduce  the  time  t, 
where  dx/dt  =  u,  as  a  parameter  and  then  x  is  a  f\mction  of  t. 
Prom  (l8),  we  have  on  a  characteristic 

y 

,2       ,,2       ^(^\_<^(y.        \  f 


{2ha) 


\^^~   1^(^^-^(^^)  +  €  J  V  Hdx 

x^ 


■■* 


v  =  v,.S%i-i-f  j 


Hdx 


where  the  integrals  are  along  the  characteristic  and  u  =  u^, 
v  =  v^,  x  =  x^;^+oo  for  t  =  t^. 

The  proof  of  lemma  2  follows  immediately  from  the  second 


line, 
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A  third  leirnna  is  needed  to  prove  the  first. 
Lemma  3:  On  a  characteristic,  the  values  of  x  where  u  changes 
from  negative  to  positive  (or  positive  to  negative)  increase 
with  increasing  time.  -     ' 

This  statement  is  equivalent  to  the  statement  that  an  ion 
has  moved  forward  each  time  it  Just  turns  around. 
Proof  of  Lemma  3: 

Suppose  the  lemma  is  not  true.  Then  for  say  x  =  x^^  we 

have,  at  some  time  t  , 
u  =  u^  =  0.  Let  u^  be  the 
value  of  u  when  the  ion 
^   first  comes  back  to  x  =  x„ 


at  time  t  =  t,. 
We  have  from  the  first  line  of  (24a), 

-—  =  €  /  V  Hdx. 

Consider  the  path  betv/een  these  two  events,  t^  <  t  <  t-j^, 

as  consisting  of  forward  (f)  and  backward  (b)  transits  (dx  >  0, 

dx  <  0)  of  an  ion. 

Then  €  /  v  Hdx  may  be  written  as  e  /  v  Hdx  -  e  /  vH|dx|. 

(f)         (b) 
Now  at  each  value  of  x  every  forward  transit  is  matched  by  a  later 

backward  transit.  By  lemma  2  the  value  of  v  is  larger  on  the  back- 
ward transit  than  the  value  of  v  on  the  corresponding  forward 
transit.   Hence  the  integral  €  J  v  Hdx  is  negative  and  therefore 
cannot  equal  u^^/P  which  is     °  positive.  Therefore  the  ion 
does  not  return  to  x  =  x^  and  the  first  half  of  the  lemma  is  true. 
The  second  half  is  proved  in  the  same  way. 
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Proof  of  Lemma  1:  Suppose  the  lemma  Is  not  true.  Then  as 


t — >  00,  X  <  X  <  00,  for  some  Ion.  From  the  second  equation  of 

^     X 
(24a)  for  t — >  00,  /  H  dx-/->  oo  since  x  <  x^.  Hence  v— >  oo. 

By  lemma  3  there  is  some  x-  and  t,  such  that  for  t  >  t^, 

we  have  x  >  x..   Suppose  u  vanishes  for  x  =  x^  >  x«  at  time  t^ 
"'  3|e|   +  k 

taken  so  large  that  v  >  — t^^ ,  k  >  0,  where  IeIj^^^^, 

mln 
H  .  J=   0  are  the  respective  maximim  and  minimum  on  x^  <  x  <  x^. 

Then  by  (24a),  on  the  first  forward  transit  after  t^ 
^4->  -   2|<t^lmax+J^5|E'max+^)^^ 

>  '^'max+^l^^-^*^ 

where  \^\  =     max   \^\  ,   k.,  is  constant. 

Hence  m  ^  0   for  x  >  x^  and  the  ion  cannot  turn  around  again. 

Therefore,  for  sufficiently  large  times  u  does  not  change  sign. 

Hence,  by  the  second  equation  of  (l8),  since  v — ^  oo,  dv/dt  = 

u  dv/dx  =  €(l-uH) — >  00,  and  thus  dv/dx  =  €(^  -  H) — >  oo  since 

u  ;^  0.  Therefore,  u^  0  through  positive  values.  But,  by  (18) 

^  duVdx  =  E  +  vH — >  00  and  therefore  u^-/-^  0.  Thus  the  hypothesis 

cannot  hold  and  the  lemma  is  proved. 
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5 .   The  Special  Case  of  Zero  Pressure. 

The  case  in  which  f^  is  a  delta  function  can  be  treated 
completely.   The  higher-than-flrst  moments  of  f°°  about  the  mean 
velocity  vanish  and  therefore  this  case  corresponds  to  zero  pres- 
sure.  We  may  also  say  that  all  the  ions  have  the  same  velocity 
for  all  values  of  x.  The  mean  velocities  u  and  v  are  easily 
shov;n  to  satisfy  (l8)  with  u  =  u  .  v  =  v  provided  u  is  positive 
along  the  particle  path.   Furthermore,  from  (23a)  and  (23b)  we 
find  explicitly  using  (l6b)  and  (12a), 

2 

u  +  ^  (1-Kn^)  =  -\ 

^   2A  2A^ 

(25) 

Av_^  +  eTl(l  +  ™)    =  0. 
From  the   first  equation  of   (24a),   and   (22) 

(26)  u^   _   1   =   -  I  ^  +  |i     /^  v^H  dx. 

—  00 

These  three  equations  yield  with  (23c)  and  (24), 

(27)  nl-   1   =\  {1   -E   -^)    -4  n^(l  +  ¥)^  . 

The  relation  between  t)  and  H  is  derived  by  eliminating  u  between 
(25)  and  (2?), 

(28)  1  +  4r  (1-H-  ^)  =  (-1  +   1  (H^-l-KTi^H^))^  4  ^^(^^  T^)  ♦ 

A^        '^  2A^  A         ^ 


The  following  solution  is  essentially  well-known  and  has  been 
established  by  many  authors,  either  using  the  particle  view- 
point as  here  or  using  a  two-fluid  approach  with  zero  pressure, 
see  for  example  [1],  [5],  [7]. 
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The  solution  curve  in  the  (H,T|)-plane  is  illustrated  in  Figure  1. 
It  exists  provided  u^  =  l  +  -L  (i_h-  ti^/2)  +  0(€  )  >  0.   It  is 
plainly  symmetric  about  ti  =  0  and  hence  the  solution  is  a  pulse, 
see  Figure  2.  There  is  no  shock.   It  is  also  clear  that  quali- 
tatively the  terms  in  K  have  no  effect. 

We  can  find  a  maximum  value,  A  „  ,  of  A  such  that  if 

'  max 

A  <  Ajj^g^,  u^  >  0  for  all  r],U   satisfying  (28).  We  may  take  €  =  0 
and  K  =  0  for  convenience.  From  (25),  for  fixed  A,  u  has  its 
minimijm  where  H  has  its  maximum  and  by  (28),  this  occurs  where 
■q  =  0.  Hence  at  the  maximiom  value  of  A,  u  =  0  where  tj  =  0. 
Therefore, 


1  +  4r  (1-H)  =0   and   1  +  -^   (l-H^)  =  0. 
A  2A'^ 


Hence  A^,^=  2. 
max 

We  shall  throughout  restrict  ourselves  to  the  case  in  which 
A  <  2.   One  can  easily  show  that  for  x  sufficiently  large  and 
negative,  E  =  -dj25/dx  <  0  and  dH/dx  >  0.   In  the  case  in  which 
f  at  -00  is  not  a  delta  function  but  a  cut-off  Maxwellian,  these 
two  conditions  will  hold  if  the  pressure  is  sufficiently  small, 
i.e.,  if  the  distribution  function  f~  is  concentrated  in  a  narrow 
region.  We  assume  this  holds. 
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6<,   The  General  Motion  of  the  Ions. 


We  proceed  to  find  solutions  of  (l8)  as  series  in  e  and 
thereby  a  formal  power  series  solution  of  (24) «  These  series 
we  write  as 


(29) 


U  =  U   +  €U,  + 
O       1 


O       1 


We  assume  E  and  H  >  0  are  given  independently  of  €,  con- 
tinuous and  bounded  for  -«  <  x  <  +»,   The  latter  condition  can 

be  verified  after  finding  an  approximate  solution.   We  find,  on 
equating  coefficients  of  powers  of  €  in  (l8)  and  integrating, 

-4) 


(30) 


2    2    2 
^o  =  ^*  -  A 


V   =  V. 


J 


■^l"o  -  J  V*  /  H  dx 

,  r  X      ^      1 

v^  =  jU   dx/u^  -  /  Hdxj 


where  u^  =  ■u(x^),  v^  =  v(x^),  j6^   =  ^{x^)   and  integrals  are  taken 

X 

along  particle  paths ^  We  write  /  dx/u^  =  t^,  the  time  that  the 
ion  takes  in  the  zero  order  approximation  to  go  from  x^  to  x. 
The  position  x^  may  be  chosen  arbitrarily  to  suit  the  particular 
application  but  we  shall  often  choose  it  to  correspond  to  a  small 
value  of  j6   and  a  large  negative  value  of  x,  small  compared  with 

l/€.  j 

The  expression  for  v^  is  not  necessarily  a  single-valued 
function  of  x.   It  assumes  a  different  value  each  time  the  ion 
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with  Initial  velocity  u^,,  v^  crosses  x  =  x^. 

The  first  order  approximation  in  €  is,  by  (50), 


2.2   2  ,jii..\    At.,    'i^  .  2e 


(31) 


yx^^   <  -  I  (^(x)-^(x^))  +  f-   vJP(x)  -  P(xJ] 

X 

v^  v^  +  f  [  /  ^ -  P(x)  +  P(x^)] 

^*  {^l-   |(j2^(x)-^(x^)) 


where  the  square  root  has  the  same  sign  as  dx  and  P(x)  -  P(xj^)  = 

X  p 

/  HdXo  The  neglected  terms  are  0(€  ).  We  must  consider  where 

X 

(31)  is  valid. 

We  note  immediately  that  the  expansion  to  zero  order  in  e 
is  not  valid  for  |x-x^|  large  of  order  l/ei  for  then,  €P(x) 
becomes  0(l), 


X 

r 

X 


Also  /  dx/>/uJ-  f(^(x)-i25(x^))  is  not  singular  where 


uj-  ^  ^  +  ^  ^^   vanishes  unless  E  ^  -d£$/dx  also  vanishes.  There- 
fore (31)  will  be  valid  for  all  ions  as  long  as  x-x^  is  bounded 
and  u^  -  2/a  {!6{x)~0{x^))   does  not  vanish  at  E  =  0. 

The  exceptional  situation  where  E  =  0  and  ^^  -  -^   ('l>(x)-<t)(x^)) 
vanishes  can  be  treated  separately  by  assuming  E  is  linear  in  x 
near  this  point.   It  can  then  be  shown  by  a  lengthy  but  straight- 
forward computation  that  the  ion  paths  either  satisfy  (31)  to 

We  shall  not  therefore  be  able  to  say  anything  about  the  motion 
of  the  particles  and  hence  about  the  fields  for  x  positive  and 
greater  than  o(l/€).  But  the  original  problem  is  not  correctly 
set  if  we  look  for  a  solution  in  the  range  -co  <  x  <  X  where  X  <  +oo 
unless  we  prescribe  on  x  =  X  the  falues  of  f^.  (and  even  of  f_)  for 
negative  values  of  u.  We  are  therefore  forced  to  impose  some  con- 
ditions for  X  large.  These  rather  weak  conditions  will  be  dis- 
cussed in  case  d)  of  this  section. 
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lowest  order  or  their  corresponding  contributions  to  the  integrals 
in  (23aa)  (23c)  tend  to  zero  as  e — >  0. 

Note  that  because  of  the  condition  on  x-x^  we  are  not  Justi- 
fied in  using  the  leading  term  alone  at  -co.  This  condition  would 
not  have  appeared  If  we  had  not  examined  the  first  order  terms 
in  €. 

On  the  other  hand  we  may  apply  (31)  when  |x-x^|  is  0(1 ) 
and  v^  is  o(l/e).  This  is  established  formally  in  the  following 
way:  We  set  E  =  E  +  ev^H,  v  =  v^  +  V.   Then  by  (l8). 


du   1  /_*  ,   ^_.\ 


dV       €    /T  tt\ 


with  u  =  u^,  V  =  0  for  x  =  x^.  These  equations  are  just  like 
(l8)  and  have  bounded  initial  conditions.   A  power  series  develop- 
ment yields 


u^  =  uj  +  I  J   E*dx  +  0(€^) 


X* 


X 


V  =  f  r      ^      p(x)  +  p(x^)  +  o{€^). 

^*  Vuf  +  I  /  E*dx 


X   *         X 


But  /  E  dx  =  /  (E  +  €V^H)dx  =  -  l6{x)   +  !6{x^)   +   €V^(P(x)  -  P(x^)) 
and  thus  we  get,  neglecting  terms  that  are  o(e). 
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26 


Vi^  =  <   -  J   (^(x)-^(x^))  +  f.  v^(P(x)-P(xJ) 


v=v^+V=v^+| 


X 

^"^  -   T(x)+T(x„) 


x*fS  -f  (^(x)-i25(xj) 


as  In  (31). 

We  shall  use  (31)  In  two  situations;  a)  where  |x-x^I  is 
o(l/€)  and  V  is  0(1)5  b)  where  |x-x^|  is  0(1 )  and  v  is  o^lA). 

We  find,  in  using  (31)  that  there  there  are  four  kinds  of 
ion  paths  to  discuss:  the  "through"  ion  paths  v.iiich  do  not 
turn  around  (u  >  0)  in  any  finite  reglor ;  the  "pseudo-trapped" 

ion  path?  which  to  lowest  order  in  e  appear  to  be  closed 

*      It        •• 
paths  ',   the  bouncing  ion  paths  which  to  lowest  order  appear 

to  be  reflected  back  to  -»  at  the  potential  peak,  and  "returning" 

ion  paths  which  appear  to  come  from  +00.   Our  objective  is  to 

find  out  what  value  of  f  should  be  associated  with  each  of 

these  paths.   (See  figure  at  the  end  of  this  section.) 

There  are  also  paths  on  which  v  is  at  least  of  order  l/e. 
They  are  not  described  by  (31).  However,  we  are  spared  finding  e. 
complete  description  of  these  paths,  as  we  shall  see,  by  the 
fact  that  f ,  =  0  there. 

a)   "Through"  ion  paths.  Prom  (31)  we  see  that  to  first 
order  in  € ,  an  ion  turns  back  (i.e.,  u  changes  sign)  only  if 

0  =  u^  =  uj  -  I  (J2$  -  4)  +  ^  v^(p(x)-P(xJ). 
To  zero  order  in  e  an  ion  follows  a  "through"  ion  path  (u  >  O) 


* 

We  recollect  from  Lemma  1  that  there  are  no  trapped  ions. 
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If  Its  initial  velocity  u^  satisfies  uf  -  f  (^  =,  "  l^i^*))  >   0 
provided  v^  «  l/e.  For  these  ions  we  note  the  important  fact 
that  V  also  remains  constant  to  lowest  order  in  €. 

b )   "Pseudo-trapped"ion  paths.  Let  us  consider  the  ion 
paths  to  zero  order  in  e  in  (u,x)-space  with  v  =  o(l/€).  Between 


u*. 


D  boundary  of  pseudo-trapped 

paths 


potential 
peak 


ps^udo- 
trapped 
path 


two  peaks  of  potential  the 
paths  appear  as  in  the  adja- 
cent sketch.  Take  any  path 
that  is  a  closed  curve  —  we 
call  the  corresponding  ions 
pseudo-trapped  ions  since  no 
ion  really  is  trapped  (Lemma 
1 ) .   We  want  to  show  that  for 


our  problem,  that  is  to  first  order  in  e,  there  are  no  ions 
(f  =  0)  on  these  paths. 

V/e  prove  this  in  several  parts.   First  we  assiome  there  is  a 
pseudo-trapped  particle,  i.e.,  such  a  particle  path  with  f_^  j^   0. 
Then  we  find  a  condition  on  the  y-component  of  velocity  for  this 
ion  when  it  starts  being  trapped.   From  this  we  find  similar  con- 
ditions when  it  first  enters  the  region  where  there  is  a  sizable 
charge  separation  firld  (E).   This  condition  in  turn  implies 
conditions  on  the  speed  of  the  ion  at  -oo.   Finally  we  notice  by  the 
cut-off  condition  that  there  are  no  ions  satisfying  this  condition 
at  -00. 

Note  that  the  ion  we  are  concerned  with  is  characterized  by 
the  conditions  x  =  x^,  u  =  u^,  v  =  v^  for  t  =  t^.   These  conditions 
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ore  Independent  of  c.   Its  earlier  position  and  velocities  are 
In  general  functions  of  t  and  of  e. 

Suppose  an  Ion  Is  near  a  point  on  one  of  these  curves  at 


(x^,u^) 


approximate 
path  P. 


some  time  with  v  =  0(1). 
That  Is,  let  u  =  u^,  v  =  v^  = 
0(1),  X  =  x^  for  t  =  t^  and 
^^  let  X,  ,  Xp  be  such  that 

^*  -  J  (i^(x)-j2$(x^))  vanishes 
for  X  =  X,  and  x  =  Xp.  Note 
that  the  curve 


u^  =  u^  -  I  (^(x)-^5(xj) 

lies  entirely  Inside  D,  see  figure. 

Vie   want  to  determine  where  the  corresponding  ion  came  from 
and  v;here  it  is  ^olng  to.   ¥e  assert  Statement  1:  This  Ion  came 
from  a  region  with  v  at  least  of  the  order  l/e. 

The  proof  Is  based  on  the  following  idea:   in  the  (x,u,v)- 

space,  the  characteristic  actually  winds  around  a  cylinder  vrhose 

cross-section  ir;  the  approximate  path  P  where  v  =  ©(l/e).   It 

differs  from  P  '.ply  for  large  v.   A  particle  moving  on  this  path 

is  In  this  region  for  a  long  time.   We  will  show  that  we  arrive 

at  a  contradiction  if  we  assume  v  is  o(l/€)  on  the  whole  path 

between  x  and  x 

o      oo 

Proof;  Suppose  v  =  o(l/e)  for  the  ion  ,   while  its  path  lies 
between  x  ,  x   up  to  the  time  t^  when  u  =  u^,  v  =  v^.   Then  on 
any  cycle  we  may  apply  (51).   If,  at  any  time  t,  prior  to  t^,  this 


ion  has  a  value  of  x^,  u,  lying  in  D,  its  path  for  t^  <  t  <  t-, 
lies  arbitrarily  close  (if  €  is  made  sufficiently  small)  to  the 
pseudo-trapped  path  passing  through  x, ,  u^.   Here  tp  is  the 
preceding  value  of  t  when  the  ion  crossed  x,  with  positive  u 
(i.e.,  one  complete  cycle).  Therefore,  to  enter  D  the  ion  must 
cross  the  boundary  D  of  D,  and  on  its  first  cycle  in  D  its  path 
lies  arbitrarily  close  to  D.   In  other  words,  on  its  first 
cycle,  by  (31),  letting  €— >  0, 

u2(x)^uj^  -  f  {S6{x)-!6{xj) 

where  u^^=  u^^(e)  is  the  value  of  u  on  D  for  x  =  x^.   Let  T(€) 

be  the  corresponding  value  of  t  with  x  =  x^  and  u(T)^^  ■^**'  ^y 

(24a),  we  have  exactly 

T 
u^(T)  =  uj  +  2e  /  V  H  dx 

where  the  integral  is  taken  (backwards  in  time)  along  the  path 
from  t  =  t^  to  t  =  T.   Hence,  combining  these  two  relations, 

T 

r 
t. 


^Mr*  ^  u^  +  2e  /  V  H  dx, 


Since  (u^,x^)  lies  inside  D  and  (u^^,x^)  lies  on  the  bound- 

*  P     2 

ary  D,  we  have  u^  <  u^^.  Thus  u^^  and  u^  differ  by  a  quantity 

that  is  0(1)  and  does  not  go  to  zero  as  e— >  0.  Therefore 

T 

e  /  V  Hdx  is  positive  and  bounded  away  from  zero  as  € — >  0.  But 
t^ 
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T  T  T      T 

e  /  V  Hdx  =  €  /  V  dP  =  e[vP]   -  e  /  P  dv. 
t^         t^  t^j    t* 

T 
Since  v  is  ©(l/e),  and  £[v]  ^  o(l),  we  conclude  that  e  /  Pdv  is 

negative  and  "boTanded  away  from  zero  as  e^   0.  But  e  J  Pdv  = 

g  r  P  ^  dx  =  €^/P(-  -  H)dx  by  (18)  or  €  /  Pdv  =  €^  /  Pdt-e^  /  F  dP 

P    "^  .2   ^  T  " 
e  /  Pdt-  ~  [P^ ]    since  dx  =  udt.  The  last  term  is  zero 

t^j. 
since  P  =  P(x)  and  x  =  x^  for  t  =  t^  and  t  =  T.  Therefore  we 

have  €  /  Pdv  =  e^  /  Pdt,  But  |€^  /  Pdt]  <  e^ !?  Imax''^"'^* '  •  ^^®^^® 
since  we  have  already  shown  that  €  /  Pdv  is  bounded  away  from  zero 

as  e — >  0,  it  follows  that  €^|T-t^|  is  bounded  away  from  zero.  By 

T 
(24a),  v(T)  =  v^  +  6(T-t^)-  €[P3  =  v^^  +  €(T-t^),  and  we  conclude 

that  €(v(t)-  v^)  is  bounded  away  from  zero  as  e — >  0.   Since 

ev^ — >  0  as  € — >  0,  by  assumption,  we  conclude  that,  as  e — >  0, 

ev(T)  is  bounded  away  from  zero  or  v(t)  y^o^l/e)  which  contradicts 

our  ass\jmption.   Statement  I  on  p.  33  is  proved. 

Furthermore,  we  shall  now  show  Statement  II:  this  ion  first 
crossed  any  plane  x  =  x^with  |x^-  x|  =  o(l/€)  at  a  time  t  =  t^with 
V  negative  if  v  i<^  ©(l/e).   We  also  derive  an  inequality  to  use 
later. 

Let  V,  't  be  the  values  of  v,  t  when  the  ion  first  crosses 
X  =  X.  We  recall  that  v^,  t^  are  the  corresponding  values  when 
the  ion  crosses  x  =  x^.  Then  on  integrating  the  second  equation 
of  motion  we  have,  with  dx  =  udt, 

v^  -  V  =  €  i  (l-uH)dt  =  €[t^-  t]-  €  I   Hdx 


or 
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11m  e(v^-  v)  =  lim  €^[t^-  t]-  lim  e^  f  Hdx. 
€^0  €-^0  €^0    -L 

Now  H  is  bounded,  and  lim  €(x^-  x)  =  0  by  ass\Amption.   Hence 

e  >0 

lim  €V  =  lim  cv^  -  lim  e^(t^-  t)  =  lim  €^(t-t^). 

Also  for  any  two  times  t^jt^  with  v  =  v^jVp,  x  =  x, ,Xp, 


(*)       lim  €(v^-  Vp)  =  lim  e^it^-   tp). 


Since  t^  >  t,   we  see  from  the  first  inequality 
lim  €V  <  0, 

€->0 

and  Statement  II  is  proved. 

By  Statement  I,  at  some  time  t, ,  with  x  =  x,  such  that 
X  <  X-  <  X   and  v  =  v,  ,  lim  ev,  ?^  0.      By  Statement  II, 
lim  €V,  <  0.   Hence  lim  ev,  <  0.   Combining  this  result  with 
equation  (*),  we  see  that  for  t  =  tp  <  t^  we  have 

(**)      lim  ev^  =  lim  ev^  -  lim  &^(t^-'t^)   <  0 

at  t  =  tp  <  t^  <  t^,  X  =  X  <  x^  where  |xp-  x^|=  o(— )  and  v  =  v'  . 

We  recall  that  this  inequality  refers  to  an  ion  which  is 
pseudo-trapped  near  x  =  x^.  We  will  next  prove  a  third  statement 
about  such  an  ion  and  show  that  this  statement  contradicts  the 
prescribed  condition  (l6a)  at  -oo.   We  pick  x  sufficiently  far 
to  the  left  so  that  for  all  e,  E  <  -a,  say,  for  x  <  x.   [In  the 
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region  to  the  left  we  cannot  use  the  approximation  (31)  since  |E| 
is  small  but  it  will  be  valid  with  appropriate  restrictions  on  v 
to  the  right , ] 

The  next  step  is  to  estimate  v  in  x  <  x^.  We  set  E  =  eS, 
X  =  xAo 

We  may  assume,  since  it  is  true  fcr  the  zero  pressure  case, 
that  £  <  0  and  d£/dx  <  0.  Then  £  changes  monotonically  from  0 
to  -a/e. 

The  equations  of  motion  become 

(  — )     "  ^  =  £^  VH,      u  §1  =  1  -  UH. 

Because  the  parameter  €  no  longer  appears  in  the  equations  of 
motion,  there  is  no  simple  expansion  procedure  for  estimating 
V.  However,  it  is  reasonable  to  treat  H  as  constant  in  this 
range  and  it  is  then  possible  to  make  various  estimates  for  the 

velocities.  These  estimates  are  given  in  Appendix  3.  We  have 

f ^ ^ 

then,  see  Lemma  A,  with  r  >  0  being  the  value  of  r  ="J(l-u)  +  v 

at  -00, 

V  >  -r 

at  each  point  on  the  ion  path''  in  the  region  where  £  <  0  and  ^  <  0» 

In  particular  this  holds  for  the  first  transit  of  the  plane 

x  =  x^o  By  (***),  it  follows  that  lim  er  >  0.  We  have  proved 

e->0 


■* 

Note  that  the  ion  in  question  is  characterized  by  the  condition 

that  it  crosses  x  =  x#  with  u  =  u*,  and  v  =  v*  where  x^,u^,v^ 

are  independent  of  e.  Everything  else  depends  on  e.  Hence  in 

particular,  r  depends  on  €. 
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statement  III;   In  the  limit  as  e -^   0,  an  Ion  which  Is  at  some 
time  pseudo-trapped  near  x  =  x^  has  at  -w  a  value  of  r  which 
becomes  Infinite  like  l/e.  But,  by  (l6a),  there  are  no  particles 
with  r  >  R,  I.e.,  f~(u,v)  =  0  there.   (Note  R  Is  fixed  Indepen- 
dent of  e.)  Hence  there  are  no  particles  on  the  pseudo-trapped 
paths.  I.e. ,  f  (u,v,x)  =  0  there. 

We  have  established,  therefore,  that  there  are  no  pseudo- 
trapped  Ions,  I.e.,  no  Ions  on  the  closed  paths  given  by 

u  +  -^  =  2  — J —  =  2  — J —  ,   V  =  constant 

with  X,  ?^  x-.   And  thus  f  =  0  there. 

c)   "Bouncing  Ion  paths".  These  are  ion  paths  which  to 
lowest  order  in  €  are  of  the  form 


u^  +  -^  (x)  =  -^  (x^)   ,  V  =  constant 


but  along  which  6{x)   ^   £5(x,  )  for  x  ?^  x,  . 

For  this  class  of  ion  paths  u  changes  sign  only  if  x  =  x^  is 
in  the  region  where  (3I)  is  valid,  that  is,  x  >  x^  say.   It  follows 
that  the  corresponding  ion  moves  with  increasing  x  from  x^  to  x, 
and  back  to  x^.   For  x  <  x^  the  magnetic  field  turns  the  ion 
around  and  it  recrosses  x  =  x^  for  some  new  value  of  u^.   At  this 
stage  the  ion  follows  a  path  which  may  or  may  not  be  of  the  above 
form.   If  it  is,  with  say  u^  +  (2/A)i25(x)  =  {2/ k)ll){x^)   the  ion  is 
again  turned  around  at  x  =  Xp  and  the  process  is  repeated.   Accord- 
ing to  Lemma  2  this  ion  acquires  an  increment  of  v  each  time  it 
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turns  aroundo  Eventually  it  passes  through  (Lemma  1),   In  this 
final  stage  we  must  have 

^^+l^>l^max 

where  't  ^     is  the  maximum  value  ot  i). 
max 

On  the  one  hand,  going  "backwards  in  time  on  the  path  of 
this  ion  we  see  that,  since  the  ion  bounces,  it  must  cross 
X  =  x^  initially  with  a  positive  value  u  less  than  >A"27a7^JJ^  . 

On  the  other  hand,  for  x  <  x^  it  can  be  reasoned  that  u 
does  not  change  on  the  average  and  that  some  particles  with 
(l-u)  +  V  =  r  at  -00  will  have  u  <  1-r  for  x  =  x^.  The  argu- 
ment is  given  in  Appendix  3  where  the  necessary  estimates  are 
developed.   Hence,  if  u  = /(2/A)0j^g^   cuts  (l-u)^  +  v^  =  r^  <  R^ 
there  will  be  some  bouncing  ions,  see  (l6a).   Of  course,  these 
bouncing  ions  are  going  to  influence  the  field  for  x  <  x^  but 
if  there  are  not  too  many  of  them  and  if  they  do  not  bounce  too 
often,  they  do  not  change  the  general  character  of  the  field  for 
X  <  x^.   This  effect  will  be  discussed  in  section  10, 

However,  their  effect  for  x  >  x^  is  of  the  greatest  signifi- 
cance as  we  shall  see. 

It  is  now  possible  to  describe  the  behavior  of  all  ions 

2    2    2 
which,  at  X  =  -oo,  have  velocities  lying  in  the  circle  (l-u)  +  v  <  R 

where  by  (l6a),  R  <  1.   They  either  pass  through  with  v  unchanged 

or  they  bounce,  possibly  at  several  places,  somewhere  between  -oo 

and  the  point  x,  where  lb   assumes  its  maximum.   It  Is  true  that 

if  they  bounce  often  enovgh,  v  may  increase  a  great  deal.  We 
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assvime     that  v  remains   0(1). 

d)   "Returning"  Ions.   These  are  ions  which  appear  to  come 
from  +00  in  the  lowest  order  approximation.  By  Lemma  1  there  are 
no  ions  actually  coming  back  from  +oo.  They  therefore  must  cor- 
respond from  our  preceding  analysis  to  ions  that  (possibly  after 
"bouncing")  went  "through". 

In  this  case  we  assume  {j>)   there  are  no  returning  ions  for 
which  u  changes  sign  between  x,  and  x.  +  o(l/€).   (This  can  be 
shown  to  be  true  if  E  =  0,  H  =  constant  for  sufficiently  large 
x,  see  figure. 


^x 


The  adjoining  sketch  indicates  the  four  kinds  of  ion  paths 

in  the  (x,u)-plane. 

Pseudo- 
Through       /trapped 


Returning 
Returning 


If  the  cut-off  radius  R  is  sufficiently  small,  no  ions  bounce. 
Hence  by  a  continuity  argvment  one  finds  that  by  increasing  R 
one  gets  first  ions  that  bounce  only  a  few  times.   For  these, 
V  remains  0(1 ) . 
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Conclusions  about  ion  paths. 

We  may  now  finally  decide  what  path  any  Ion  follows  from 
-00 o   An  ion  may  pass  through  the  whole  region  to  large  positive 
X  without  turning  around.   It  if  turns  around  it  bounces  "back  to 
large  negative  x  but  eventually  passes  through.   On  its  whole 
path  V  is  o(l/€).   An  ion  never  gets  on  a  pseudo-trapped  path. 


7,  Estimating  ev  and  Determining  the  Potential  t>. 

Now  that  we  know  the  behavior  of  the  ions  we  can  apply  our 
knowledge  to  find  how  the  fields  behave. 

On  through  ion  paths,  v  is  approximately  constant  by  (31) 
and  hence  v  is  0(1).   On  bouncing  ion  paths,  v  does  not  change 
by  an  order  of  magnitude.  There  are  no  pseudo-trapped  ions  and 
by  the  assvimption  for  returning  ions  v  is  o(l/€).   Hence  by  (11), 
V  =  //  f.v  du  dv  /  //  f ,du  dv,  we  see  that  €v,- — >  0  as  e — >  0, 
Therefore,  from  (23c), 

(3?)      Aj6  =  constant  +  H  +  tj  /2, 

We  also  note  that  11m  J  =  0  and  lira  ct]  =  0.   Hence  (23b) 

€^0  €^0 

yields  no  Information  in  the  lowest  order  case. 
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8.   Boltzmann  Distribution  for  Ions. 

Having  discussed  In  Section  6  the  general  motion  of  the  ions 
we  are  now  able  to  determine  f  to  zero  order  in  €  for 

-00  <  X  <  o(l/€).  Then  in  the  next  section,  we  shall  discuss  the 

2    2  V 
fields.  By  (31),  u  +  j  0  =  constant  on  ion  paths  where  €v^  0. 

Hence  f  (x,u,v)  is  constant  on  u  +  -^  0(x)  =  constant.   Hence 
f  is  even  in  u  for  those  values  of  u  and  x  for  which 
u  +  y  J^(x)  <  J  !6{x^).      On  X  =  x^,  f  is  determined  on  all  particle 
paths,  which  meet  x  =  x^  for  the  first  time,  by  its  value  at  -oo 
on  the  corresponding  particle  path.  Because  of  bouncing  and  return- 
ing, certain  relations  among  the  values  of  f  for  some  positive  u 
hold,  which  are  given  roughly  by  taking  H  =  1,  t)  =  0  for  x  <  x^, 
if  x^  is  small  enough,  but  we  do  not  need  to  know  exactly  what  they 
are.  We  write  for  x  >  x^ 

f^(x,u,v)  =  f^{u,v,6)   =  F^{   u^+   ^  ,  v)  for  u  >  o 


(32a) 


oru2+^<|^(x,) 


Ct^i   u^+  ^  ,  v)  for  u  <  0  and  u^+  ^  >  |  !6{x^) 


and  these  expressions  will  be  valid  as  long  as  ^   is  monotonically 

increasing.   Where  s6   reaches  its  first  maximx;im  at  x  =  x^  say, 

the  ion  paths  u  +  x  ^  =  constant  appear  as  in  the  figure.   We  will 


^x 


x=x, 
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show  below  that  f  =  0  In  the  region  D  to  the  right.   In  the  region 
A  to  the  left  f  Is  determined  by  the  value  of  f_^  at  x  =  x^  and 
may  or  may  not  be  zero. 

To  show  f  =  ^'  In  D  consider  the  two  possible  cases:  either 
ii   has  its  first  minlmiam  for  x  >  x,  at  x  =  X  for 

(1)  X  larger  than  o(l/e)  or 

(2)  X  =  x^^  =  o(l/€). 

In  the  first  case  we  have  the  situation  for  which  the  third 
assumption  about  returning  ions  (case  d)  was  designed.  There 
are  no  Ions  in  D.   Hence  f^  =  0  in  D. 

In  the  second  case  let  x,p  and  x,,  be  two  values  of  x 
straddling  x,^  and  such  that  ^^{x^g)  =  ^^^t^^  ^Wa   E  =  -  dj2$/dx 

positive  for  x,p  <  x  <  x,^ 
and  negative  for  x,,  <  x  <  x^^. 
Then  f ^{yi^^.yx ,v)   =  f_^(x^^,u,v) 
since  i{y.^^)   =  j2$(x^^).   From 
(23aa)  and  {23a)  we  see  that 
I  and  hence  H^(I-Kti^)  have  the 
same  value  at  x-,p  and  at  x,^. 

Also  since  ^^^(x^^)  =  ^(x^^,),  by  ((32),  p.  4l),  H  +  ti^/2  has 

2 
the  same  value  at  the  two  points.   Hence  H  and  .-i  are  the  same 


for  all  such  pairs  of  points.   It  follows  furthermore  from 

X         , 
j6(x)-j?5(x^^)=  -/  Edx  =  /  J  Hdx  that  Ix^^"  ^n '  =  l^^^l^"  ^11 ' ' 

Hence  the   curves  in  the  x,u-plane   on  which  u2  +  2j?$/a   is   constant 


are  symmetric  about  x,,  and  therefore  D,  which  is  bounded  by  the 

curves  u  +  2tf k  =  ^  i  ^      is  closed  at  both  ends  since  it  is  clost 

A  max  ^ — 

at  one.  Therefore,  inside  D  the  ions  are  "pseudo-trapped"  and 
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hence  f  =  0  In  D. 

We  are  now  finally  able  to  associate  with  each  character- 
istic to  lowest  order  in  €  a  value  of  f  .   If  the  characteristic 
comes  from  x  =  x^,  x^  large  and  negative,  the  value  of  f  on  the 
characteristic  is  the  value  at  x  =  x^,  and  now  we  have  shown  that 
if  it  does  not  cut  x  =  x^,  f  =  0. 


9.  Behavior  of  the  Fields. 

We  now  proceed  to  establish  certain  inequalities  through 
which  we  will  be  able  to  draw  qualitative  conclusions  about 
the  fields. 

We  consider  I  as  a  function  of  x,  see  (235aa)  and  (3?  ), 

^     rr    2 

I  =  //  f^u'^du  dv 


(33)    <        =  //  f^  (u,v,j?5)u^du  dv 


Vv 


for  x^  <  X  <  x. 


=  Ii(^) 

where  x,  is  the  first  maximum  of  ^,   see  (32a).   Similarly 
I  =  //  f_^u^du  dv  =  //  fg  (u,v,j2$)u^du  dv 


(3M 


^  =  l2(^) 


for  x^  <  X  <  Xg 


where  Xp  is  the  second  maximum  of  !6.      We  note  that  f,  (u,v,j^)  > 
fp(u,v,j?5)  since  they  are  either  equal  or  fp  =  0.  For  through 
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ions,  f^  =  fg.   Hence  1^(6)   >  l^ii^)  ■>      In  particular  also  f^  =  f^ 
for  l6  =  s6^   and  thus  I^^CJ^]^)  =  IgC^l)-   ^^  follows,  for  x  >  x^, 
that  we  must  distinguish  between  two  cases  according  to  the  value 
of  the  cut-off  radius  R,  see  (l6a)o 

Case  i)  There  are  no  bouncing  ions  and 


Case  ii)  There  are  bouncing  ions  and 
I^(j^)  >  I^{6). 

We  can  now  discuss  the  fields » 

In  case  i)  we  see  from  (23a),  (S^aa),  and  (32)  that  each 
pair  of  values  of  x,  one  greater  and  one  less  than  x  with  the 
same  value  of  j6   also  have  the  same  value  of  H  and  t)  and  hence, 
by  (32)  and  dj^  =  -  Edx  -   t|H  dx/A^  the  same  value  of  [x  -  x,  [. 
Hence  the  solution  is  synurietric  about  x  =  x-,  and  is  a  pulse.   In 
the  (ti,H) -plane  the  solution  is  similar  to  that  given  in  Figure  1 
(Section  1),   There  is  no  shock.   This  situation  can  be  achieved 
by  taking  the  cut-off  radius  R  small  enough  (see  (l6a))  since 
the  solution  will  be  as  close  as  we  like  to  the  delta -function 
solution  discussed  in  Section  4. 

Case  ii)o  In  this  case  we  see^,  using  (32),  (23a)  and  (23aa) 
that  at  two  values  of  x  straddling  x,  ,  say  Xp,  <  x,  and  x^,  >  x,  , 
corresponding  to  the  same  value  of  H,  we  have 

l^ilSix^^))-  KR^!6{x^-^)/A   =  l^{l6{x^^))-   KH^j2$(x^^)/a, 

Since  along  a  path   u    decreases  with  l6   increasing,  see  (31)^ 
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l-^{^)   and  IgC,?^)  are  decreasing  functions  of  ^.      The  same  is 
true  of  l^{l6)-   KH^  !6/a   and  I^{j6)-   KH^  iS/a   for  fixed  H.   Hence 
since  l-^{s6)   >  1^(6),   !6{x^-^)   >  !6{x^^) .      In  turn,  by  (32), 
•n  (x^-j^)  >  r\   (x^^).   Therefore,  in  the  H,Ti-plane  the  solution 
curve,  for  x  >  x, ,  lies  inside  the  reflection  of  the  solution 
for  X  <  X,  about  t]  =  0,  see  figure. 

x<x. 


The  arrow  indicates  increasing  x. 
In  particular,  we  may  find  a  value  of  H  -  1  so  that  the 

corresponding  r;  is  sufficiently  small  with  x  <  x^,  and  t|  for 
x  >  X,  vanishes.   At  this  point,  x  =  x^,  l6   has  a  minim\;im.  By 
symmetry  and  arguments  similar  to  those  of  Section  8,  it  can 
be  easily  shown  that  for  x  >  x,  ,  I  =  Ip(j2$)  and  hence  the  solu- 
tion remains  symmetric  about  rj  =  0  for  x  >  x,  . 
Character  of  fields. 

If  we  look  at  the  whole  field  over  x^  <  x  <  +oo  we  see 
that  the  magnetic  fie.Td  H  increases  and  reaches  a  maximum 
value.   It  then  oscillates  between  this  maximum  and  some  mini- 
mum that  is  greater  than  1. 

It  is  planned  to  find  from  computations  with  distribu- 
tions that  are  cut-off  Maxwellians  at  -co  the  solutions  for 
various  cut-off  radii  R,  gas-dynamic  Mach  numbers  M  at  -co.  and 
Alfven-Mach  numbers  A.   A  simple  estimate  shows  that  the  minimum 
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value  of  VLZ   -  2s6/a   with  f .  ?^  0  decreases  linearly.   Hence  case  11) 
will  always  occur  for  fixed  A  and  M  if  R  Is  Increased. 

In  these  computations  It  Is  Impossible  to  handle  the  charge 
separation  field  for  x  <  x^o   Therefore,  we  simply  take  H  = 
constant  and  E  =  0  (or  t)  =  O)  for  x  <  x^.  The  lowest  order  theory 
for  X  >  x^  corresponds  to  taking  H  =  0  for  x  >  x^  for  the  ions. 

A  third  case  is  perhaps  worth  mentioning  although  it  is 
Incomplete. 

Case  ill)  With  the  above  additional  approximation  for 
large  negative  x  we  find  that  an  ion  which  arrives  with  u  =  0, 
V  >  0  at  X  =  x^  remains  there.   It  cannot  move  forward,  see  (pi) 
with  €  =  0,  since  y-2j5/A  is  imaginary.   It  cannot  move  back- 
wards  into  the  constant  H  field  since  u  =  /  €vHdx  becomes 
negative.  Hence  the  approximation  cannot  be  valid.   However, 
this  bad  approximation  is  still  illtjminatlng.   We  choose  a  cut- 
off radius  so  that  no  particle  is  stuck  at  x  =  x^  in  this 
approximation  but  t\arns  around  a  great  number  of  times  before 
excaping.   This  means  it  acquires  a  very  large  value  of  v.  There- 
fore, the  distribution  function  for  x  >  x^  has  a  hump  where  v  is 
very  large.   In  other  words,  large  distortions  of  the  Maxwelllan 
for  large  velocities  are  produced  by  a  niomber  of  ions  which 
started  with  small  velocities.  However,  we  have  found  numerically 
that  if  we  let  R— ^  1,  not  only  does  the  distortion  become  infinite, 
but  also  the  pressure  and  current,  so  that  the  conditions  for  using 
a  lowest  order  theory  are  violated.   We,  therefore,  need  a  special 
theory  for  these  ions  and  this  has  not  been  worked  out  as  yet. 
Meanwhile  the  method  should  be  applied  only  to  cases  where  the 
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ions  do  not  turn  around  too  often. 

However,  the  singularities  in  the  pressure  and  current 

Indicate  that  some  of  the  ions  which  at  -oo  have  R^l  play  in 

the  full  problem,  i.e.,  where  no  approximations  are  made  for 

X  <  x^,  a  special  role  which  should  be  analyzed  but  is  probably 
qualitatively  as  described. 


10.   Actual  Behavior  of  Fields  as  x  — >  -oo. 


Certain  properties  of  the  fields  for  large  negative  x  have 
been  assumed.   We  would  like  to  sketch  here  an  asymptotic 
development  for  this  region  which  seems  to  be  consistent  with 
our  present  knowledge. 

We  want  to  combine  two  limiting  processes,  letting  x — >  -oo 
and  € — >  0,  and  we  must  do  this  in  the  right  order.   If  we  let 
€ — >  0  we  see  that  the  point  at  -co  corresponds  to  l/u  =  H  where 
u  =  (//f  du  dv)~  ,  see  '12)  and  (12a).   If  we  have  bouncing 
ions  and  evaluate  u  to  lowest  order  in  e  for  x„ — >  -co  we  find 
that  "u  (-oo)  <  1  and  hence  H(-oo)  >  1,  which  contradicts  the  pre- 
scribed conditions.   The  correct  process  for  finding  the 
behavior  at  -oo  is  to  let  x — >  -oo  and  then  e — >  0.   We  may  regard 
our  lowest  order  solution  as  being  a  good  approximation  only  as 
long  as  |-x|  «  l/e  and  E  »  €,  see  (31).   This  situation  is 
completely  analagous  to  the  usual  boundary  layer  theory. 

One  cannot  say  very  much  about  the  behavior  for  |-x|-^l/€, 
but  some  qualitative  remarks  and  justifications  may  be  made.   We. 
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return  to  equation  (21)  and  for  convenience  let  us  set  K  =  0. 
We  stretch  certain  of  the  variables »  We  set  tj  =  €Y,  and 
ex  =  X.  Then  we  have 


dH  .  y(H.c^  ^)      '"+ 


dx 


dX' 


u 


0  =  £-  -  H  +  £ 
u. 


2  dY 
c[X 


Sf, 


Au 


+ 


c^f, 


^+  (.  ^+v)H^+  (1  -UH)^  = 


0. 


Next  we  let  e- 


0  and  find 


(35) 


1  ^  =  Y  +  A  V 
H  HX    ^  +  «v_j_ 


-    1 

"+  =  H 


M,  y  ^f .  c^f , 


=  0. 


Now  if  no  Ions  turn  around  v  is  small  (it  is  zero  at  -oo 
and  small  at  x  ==  x^)  and  can  be  neglected.   This  is  confirmed 
in  the  zero  pressure  example  where  v  =  0(€  ).   But  if  there 
are  ions  coming  back  into  the  region  where  (35)  is  valid,  from 
higher  values  of  x,  v  will  be  positive  since,  1)  by  Lemma  2,   v 
Increases  at  x  with  time  and,  2)  the  contribution  to  v  from  the 
ions  on  their  original  passage  to  this  value  of  x  is  small. 

Let  us  suppose  we  have  a  zero  order  solution  as  described 
In  Section  9  in  which  H — ->  H  as  x— ->  -oo.   In  the  r|,  (H-1) -plane 
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and  Y, (H-l) -plane  we  have  for  z  ^  0   the  solid  curves,  while 

If  we  let  € >  0  we  will  obtain  the  dotted  curves  which  look 

quite  different  in  the  two  figures. 

^  increasing  x 


This  "behavior  checks  qualitatively  with  both  the  zero 
order  solution  and  (35)-   Suppose  we  are  approaching  H  =  H 

on  the  solid  curve,  e  ^  0,  with  decreasing  x.   We  are  within 

1  dH 
€  of  the  dotted  curve  in  the  (Ti,H-l)-plane  and  ^  ^  — >  0 

until  T|  is  0(€).   Then  we  must  use  (35)  and  we  see  that  Y  is 


not  going  to  zero  like  dH/dx  but  remains  larger  because  Av  >  0. 
On  the  other  hand,  the  density  is  decreasing  because  gradually 
we  reach  a  region  where  no  more  ions  are  turning  around.   There- 
fore, l/u  is  decreasing,  or,  by  (35),  H  is  decreasing. 
Eventually,  we  reach  a  region  which  all  the  ions  pass  only  once, 
the  magnetic  field  is  nearly  constant,  and  the  charge  separation 
field  is  almost  zero.   The  fields  then  approach  their  constant 
state  at  -oo  in  accordance  with  the  behavior  described  by  Calkin 
[3]. 
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APPENDIX  1  -  Unit  of  Length 

We  shall  relate  the  unit  of  length  introduced  in  (8)  with 
the  distances  the  ions  and  electrons  move  forward  in  a  single 
period. 

The  equations  of  motion  of  an  ion  are,  in  m.k.s.  units, 


d  x'   ^  du'    ^f~.         ,,„iij   N 


dt 

,2 


^-Z;  =  m^  ^  =  e(E,,,  +nu'H,,) 
dt 


"^+  TT?  =  '"+  d^ 


With  E  .  =  0,  E  ,  and  H  ,  constant  (and  evaluated  at  -co),  we 
x'      y'      z' 


find 


u'-Ey,/^H^,  =  (u--Ey,/M.H^,)  cos(^^iH^.t)-v^  sinCj^nH^.t) 


v-  =  (u^-  Ey,/nH^,)  sin(j|-|iH^,t)  +  v^  cos(jj^  ^H^.t ) 

where  u',v'  are  the  values  of  velocity  at  t  =  0.   The  time  T 
o  o  ■•■ 

for  the  ions  to  pass  through  a  complete  period  is  2Tm_^/e\iU^, . 
The  distance  i  an  ion  moves  in  T  is 

7+       2TTm   E  , 

f    u'dt  =  __t  _y_  . 

o  "^  z' 


Similarly  the  distance  i_  an  electron  moves  forward  in  T_  = 
(27nn_/e)(Ey,/nH^,)  is 
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o     E„,  2iTs/^       E    ,      /m     n 


=  2Tr  —     — =  27rA  — 

lAl  H_, 


by  (8)  and  (9)«   Since  in  all  cases  considered  1  <  A  <  2  we 
see  that  in  the  dimensionless  system  (8),  (9)  length  is 
measured  in  units  of  I,     This  Is  then  the  characteristic  length 
of  the  motion  and  the  width  of  the  transition  zones  or  the 
periods  of  the  oscillation  in  case  ii )  will  "be  some  multiple 
of  I   and  if  there  are  enough  bouncing  ions  it  will  be  of  the 
same  order  as  i. 
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APPENDIX  2  -  Asymptotic  Behavior  of  Electrons 

To  determine  asymptotic  expressions  for  the  mean  electron 
velocities  we  transform  the  Boltzmann  equation  for  the  electrons 
using  variables  given  In  (l8): 


=  u  -  u 


4  1 

=    U    ■ 

1 

^ 

H-S 

Tl' 

=    V    ■ 

■?.    . 

f  (x,u,v)  =  F  (x,^',Ti'). 


Then  the  equation  for  P  Is: 

A (u  +  e '  )  ^  -   [ A(ir  +  I .  )g  +  V^]  ^  +  ^  ^  =  0. 
setting  P  =  H^/^  u-l/\   Q  =  u'^^^  H'^/^ 

C  =  p(x)fi  cos  e 

T)'    =  Q(x)fi  sin  0 

P   (x,e',n')   =  P(x,fi,a) 


we  find  with  the  notation  u  '^—  =  t  ,  -^  =  t^ 

oX         cX     X 


"  ^3r  -  ^    2  ^P  +  Q^  ^  -  2  ^P    Q^  ^  ^°^  2^ 


.  I  (J  -  |)  II  sm  20 
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+  T-  W  -  Q-)  ^  ^^"  ^ 
_^  PQ  /  ^x  _  \\  BF  ^.„  ,f, 


-  -   COS  0  ^F   sin  e   ?.F 
^1      P   ^  "  Pfi   ^ 


r^F        ^   -^  ^ 


1  fP   Qn  ^  sin  26 


^P    4  ^^-^    1  .2  SF 


J(nH^-i^)=|p 


u-     -    ^ 


P    p2   HQ 

since  —  =  P  =  -p-  • 

Adding  all  the  above  expressions  together  and  noting  that 
du/dx  =  -(^  +  |)  we  find  the  differential  equation  for  F(x,n,©) 


is 


CJF  .  Va  g  cos  9    -tPij2(|  ^  +  1  ^)  .  ^  ]  g  COS  9 
P     0      ^     /^  '    ^^ 

,      r?Q     f     X  ^Xn   ,  U  t  SF  „^.„  q   .  f^Q  oF  ■   ._  ofl 
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^   Pfi  ,  X   "^x  N  SF  si^  -.0  C   p2  c5F  . 


<\;      /V 


We  may  now  find  an  asymptotic  expansion  for  F  =;  F  +  €F,  -t-  , 
by  substitution  and  integration  with  respect  to  0: 

Fq  =  T(x,J7) 

^1  =  ^  ?(x,fi)0  +  Pfi  II  sin  0  -  [PQ^   (|  -21  ^  1  _2i  ) 

Since  F,   must  be  a   single  valued  function  of  0,  7=  a  v^  =  0, 
We  may  choose   J,    =   0.      The  mean  velocities  are   found   using    (11), 
(12)   and  du  dv  =  du'dv'    =  PQfidfid0, 

u  du  dv  =  PQfi(u  +  Pfi  cos   e)dnd0 


V  du  dv  =  PQfi(Qfi  sin  0  +  ^)dnd0, 


We   find 


V.    =  ;^+   0(e)    =    -  I  +   0(6) 


"  +  d3r  ^TT^ 


as  required  in  (20) 
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APPENDIX  3  -  Estimates 

In  this  appendix  we  shall  get  some  bounds  on  the  velocities 
of  an  Ion  In  a  constant  magnetic  field  but  a  decreasing  charge 

separation  field.   In  equations   (l8)  we  take  (Al)  H  =  1,  E  <  0, 

dF 

^  <  0,  E{-oo)  =  0.   Equations  (l8)  become 


(Al) 


^=u^  =  -i(^+v) 
dt   "  dx   A  ^  ^  ^^ 


dv     dv   e  /n   V 

dt  =  ^  d3r  =  A  (^-") 


These  equations  we  regard  as  equations  of  motion  for  a  particle 
path  in  u,v,x-space  with  t  as  time.   Then  from  (l8)  we  find 

(A2)      ^t{u-l)2+(v  +  |)']  =  2(v^|)  .  Ig 

or  integrating  along  an  ion  path  to  a  point  P  from  some  point 
P  with  X  =  -00  we  have,  using  Al, 

00 

(A3)  {u-lf  +   (v  +  |)^  =  r^  +  2  /(v  +  |)  1  g  dx 

00 

where 

(kk)  r^  =   [{u-lf  +  v^] 

CO 


The  Integral  can  be  estimated  first  for  certain  P,  namely 
those  which  correspond  to  a  first  transit.   We  say  P  belongs  to 
a  first  transit  if  for  t  <  t(P),  x  <  x(P).   The  result  is 


In  general  it  is  sufficient  to  take  H  constant,  dE/dx  of  fixed 
sign. 
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Lemma  Al ;   If  P  is  a  first  transit,  v(P)  >  -r. 
Proof;  /  (v  +  — )  ^rr  dx  may  be  written  as  the  sums  of  the 
separate"*  integrals  along  the  forward  transits  and  backward 
transits  of  the  Ions,  i.e.. 


p 


(A4) 


/(v  .  f  )f  dx  =  ^  /  (V  .  f  )g  dx  .  H  ;  (V  .  |)||  dx 
P      ^  ^     forward     ^  ^^    backward     ^  ^ 

00 


=  ^  /  (V .  |)g  idxi  -  y:  /  (V  +  |)f  itoi 

forward  ^  ^^  baclcward  ^   ^^ 


To  each  backward  transit  there  corresponds  an  Immediately 
subsequent  forward  transit.  The  first  part  of  this  forward 
transit  has  points  which  correspond  to  points  on  the  backward 
transit  and  this  is  followed  (see  lemma  3)   by  points  which  are 


kT^ Points  of  first 
transit 


ackward 
transit 


of  first  transit,  see  figure.  Thus  the  integral  can  be  rewritten 
as 

Av  .  |)S  dx  =  5  [  /(v  .  |)g  Idxl  -  /(v  .  |)§  Idxl  ] 


(A5) 


1   Pi  B^ 


+  /   (v+|)^  Idxl 


S"(P)  "    ^^ 
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where  S(P)  Is  the  set  of  points  of  first  transit  prior  to  P, 
B.  is  the  set  of  points  of  the  1   backward  transit  and 
F  is  the  set  of  points  of  the  immediately  subsequent  forward 
transit  that  are  not  in  S.   Now, 


/(v.|)g|dx|  -/(v.|)gldx|  <0 
^i  ^i     . 


because  the  left  hand  side  is,  in  obvious  notation. 


/  (vp  ~  "^B  ^  dSc  I^"'^'  over  the  x  interval  determined  by  B.  , 

dE 
But  by  (Al),  :5—  <  0,  and  by  lemma  2,    Section  4,  v„  >  v^  . 
ux  if .     D. 

Hence  the  integral  is  negative.   Therefore 

(AS)       /("  ^  |)S  ^^  <  /(v  ^  f)g  i<^xi  <  /  V  g  idxu  s!m 

00 

by  (Al).   Substituting  (A6)  in  (A3)  we  have 

(A7)      (u-l)24.  (v  +|)^]   <  r^  +1  /  vg  |dx|  +4 


or 


a 


€    p         €  3(p)   dx  ^d 


Let  V  be  the  minimum  value  of  v  on  S  and  let  it  be  attained 

dE 


t   P  =  P^^.      Then  since  E  <  0,  g  <   0,    see    (Al), 


(V.|)'<r2.|     /        v§|dx|    .4<r2^|vE.4 

^  ^    A(P^,)       ^  €  ^  €'^ 


_    C 


7   - 


or 

(a8)      V^  <  r^ 

and  the  lemma  is  proved. 

The  second  estimate  concerns  estimates  for  points  that 
are  not  points  of  first  transit. 
Lemma  A2;  For  all  points 

o  T7  2        |eL  2   p     2|EL,pv 

(u-l)^(v  4-  f)  <  (r  +  yP)  +  |(r  +  -iyiPl)(|E|    -|E|  ), 
€   -        E        e  E  ^(p)     p 

where  ^(P)  is  the  maxlmijm  value  of  x  on  the  characteristic  up  to  P. 

In  (A5)  we  note  that  the  end  point  of  the  integration  may  be 
in  P.  or  B. .  We  find  by  similar  arguments  to  those  used  above, 

(A9)      /  (V  .  f)§  dx  <  /  (V  +  |)g  Idxl  -  /  (V  .  |)g  idxl 

CO  ^ 

where  L  is  that  part  of  the  last  backward  transit  before  P  which 
is  not  covered  by  a  subsequent  forward  transit.   Now,  by  (a8),(A1), 

(AlO)     /  (v  +|)g  Idxl  <  r|E|  +1^  ] 

S      ^  ^"^  ^^  x=^(P) 

where  ^  is  the  maximum  value  of  x  on  the  path  of  the  ion  up  to 

P.   On  the  other  hand,  by  lemma  3,  the  value  of  x  at  the  right  extreme 

point  of  L  will  be  ^(P)  and  at  the  left  extreme  will  be  x  =  x(P). 

Hence 

(All)     -  /  (V  .  |)f  Idxl  =  -  /  V  g  Idxl  .  |!  ]^.  |!  ] 


Thus   combining   (A9) , (AlO) , (All) ,   we  have 


(A12) 


/(v.|)Sdx<r|E|.|!]^-/vf   Idx 

00 


Substituting   (A12)   into   (A3)  we  find 


(A13)  (v  +  |)   <   (u-l)2+(v  +  |)   <   (r  +  M)      -  I  /  V  g   |dx| 

L 


Let  V  <  V  on  L  and  let  v  =  V  for  some  point  P, .   Then,  by 


(A13) 


2    2   2r  E    2 
—       €     e 


2r|E| 
"e" 


2   '^■'l^l?    2 
+  -  IE  V  <  r  +  — ^ — ^  +  4  |EL  •  V. 


€  '^'^ 


Hence 


or 


V  - 


< 


(r  +  -z-^-) 


(Al4)     V  <  r  + 


2|E| 


on  L. 


Substituting  (Al4)  in  the  second  inequality  of  (A13)  we  find  for 
all  points 

?      !  K I   ?         ?  I  p  I 
(A15)    (u-l)2+(v  +  |)  <  (r  +^)   +1  (r  +— ^)(|El^(p)-|E|p) 

and  the  lemma  is  proved, 

2.   Behavior  of  ions  that  do  not  turn  around. 

Here  we  wish  to  give  reasons  why,  when  an  ion  passes  through 
the  region  where  E  is  small,  the  value  of  u  does  not  change  by 
an  order  of  magnitude.  We  are  concerned  here  with  ions  that  are 
not  turning  around  in  this  region.   (They  might  turn  around 
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where  E  gets  larger. ) 

The  inequalities  of  Lemma  A2  do  not  help.   We  need  exact 
information  about  the  trajectories  which  is  impossible  to  get 
in  general.  We  therefore  consider  the  case  in  which  E  changes 
linearly  from  0  to  -a,  say,  over  a  very  large  distance  so  that 
we  may  consider  H.   =  jdE/c  dx  |  small  compared  with  1.   Again  H 
is  taken  to  be  unity  without  loss  of  generality.   We  set 
E  =  -a  -cilx  over  -a/ i   <.  x  <  0,  where  i  >  0.   Then  the  equations 
of  motion  are,  by  (l8), 

where  dx/dt  =  u. 

Differentiating  with  respect  to  t  we  have 

^  =  -  €iu  +  €(l-u) 


dt 
and  thus  we  find 


u  -  ^  =  [u  -  ^  ]  ,  cos  €  /in  (t-t(^)) 


T+I  -  ^"  -  iTT  Jg/^  "^^  ^  '^^"^   \---\-jr^ 


+  [v]    sin  €  \/l+I  (t-t(-a/i)). 
Setting  u  =  1+r  cos  |,  v  =  r  sin  ^  at  -a/ I   we  see  that 


(A15)     u  =  1  +  —J-   [cos  e  l/m  (t-t(-a/^))-l] 

+  r[cos  e  l/l+I  (t-t(-a/i))-|]. 
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The  second  term  is  always  negative. 

Hence  the  effect  of  the  charge  separation  field  is  to  change 
the  frequency  of  the  motion  and  to  slow  the  particle  down.  That 
is,  u  is  slightly  decreased,  since  I   is  small,  from  its  range  of 
values  without  the  charge  separation  field. 

Furthermore,  from  (AI5),  at  x  =  0  we  see  that,  if  u-1  is 

/        2    2 
close  to  r  at  X  =  0  the  ion  had^  at  x  -  -a/i,  (1-u)  +  v  close 

to  r  . 
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